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Nowadays, the Cosmic Microwave Background (CMB) anisotropy studies accurately 
determine the baryon fraction Ub, showing an overall and striking agreement with previous 
determinations of ujb obtained from Big Bang Nucleosynthesis (BBN). However, a deeper 
comparison of BBN predictions with the determinations of the primordial light nuclide 
abundances shows some tension, motivating an effort to further improve the accuracy of 
theoretical predictions, as well as to better evaluate systematics in both observations and 
nuclear reactions measurements. We present some results of an important step towards 
an increasing precision of BBN predictions, namely an updated and critical review of the 
nuclear network, and a new protocol to perform the nuclear data regression. 

1. Introduction 

In the framework of the "Cosmological Concordance Model" , BBN probes the earliest 
times, but once ujb is fixed, as now allowed by accurate CMB anisotropies analysis [ lj, it is 
completely ruled by Standard Physics, at least in its minimal formulation, thus providing 
useful insights on a wide range of astrophysical and cosmological issues. 

As well known, the nuclide abundance predictions are mainly affected by the nuclear 
reaction uncertainties. The suggested range of u>b and the low-energy nuclear data taken 
in the last decade justify a revision of the BBN network reliability, with respect to the one 
performed in the seminal paper by Smith et al. here we summarize our techniques and 
some results. For a detailed discussion of the issues presented here and of other refinements 
(like a better treatment of the neutrino sector and of plasma and QED effects) see 



2. From the data to the rates 

The reaction rates are obtained as thermal averages of the relevant astrophysical factors 
S{E) 1 , so the first and most critical step is how to combine the data {ik} of several 
experiments {k}, each one affected in general by statistical, <7j fc , and normalization errors, 
et, in a meta-analysis where both the magnitude and the energy behavior of the S(E) 
have to be deduced. Usually, different experiments disagree within the quoted errors, 
suggesting some systematic discrepancy, the bulk of which can be attributed to different 

^dn what follows for simplicity we will refer to the S-factor, clearly defined only for reactions induced by 
charged particles. Similar relations also hold for neutron induced reactions. 
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normalizations. There is obviously no unique and unambiguous way to deal with such 
discrepancies and several methods appeared in the literature. In a recent one described 
in [H], both the errors are included in the covariance matrix entering the expression of the 
X 2 function to be minimized. This approach clearly takes into account the correlations 
between data of the same experiment, but as addressed in [E], a bias is thus introduced, 
leading to a systematic underestimate of the fitted functions, that typically pass below the 
majority of the data, and possibly affecting the nuclide predictions and errors estimates. 
In our method, that generalizes the approach suggested in [E], the x 2 is calculated as 

v 2/_ , , \ _ v 2 .2 _ \- (Sth(E ik ,ai) - S ik u k ) 2 ^ (u k - l) 2 

X [ai, UJ k ) - Xstat + Xnorm = 2^ m h }^ 2 ' 

where Si k is the data point at the energy E ik and Sth(E) is the theoretical value of the 
astrophysical factor, depending on some parameters ai to be determined together with 
the renormalizing factors uo k by standard minimization procedures. The covariance matrix 
was built by considering the statistical errors only 2 , while different renormalization fac- 
tors were allowed for each data set; the introduction of Xnorm disfavors renormalizations 
greater than the estimated e&. The best-fit curves thus produced pass through the data, 
closer to the best determinations, as one would expect from an unbiased estimator (see 
Figs. 1 1121) . In such an approach, the bulk of the systematic uncertainty has been taken 
into account and the residual discrepancies can be considered as due to some uniden- 
tified/underestimated source of error in one or several experiments. Then, we simply 
inflated the calculated error by a scale factor -y/x 2 , as prescribed by the Particle Data 
Book [ 6J. Using this prescription, the overall normalization uncertainty cannot be clearly 
worse than the one determined by the most accurate experiment: this error was quadrati- 
cally added to the previous one. Finally, a typical renormalization factor e was calculated, 
to give an idea of the current disagreement on the scale of S(E) among several experiments 

where w k = (xt/Nk)' 1 , N k is the number of data and xl the contribution to the x 2 °f 
the k-th data set, thus assigning a larger weight to the experiments closer to the fitted S. 

The rate R was obtained by numerical integration of S(E) convolved with the Boltz- 
mann/Gamow Kernel K(E, T) ~ e - E / T e ~V E G/ E 

R= / dE K(E, T) S(E, a) (2.3) 
Jo 

and its error SR trough the standard error propagation as 

5R 2 = r d E'K(E',T) r dEK(E,T)Y, d -^ ) 
Jo Jo ft da L 



da r 



cov(ai,a m ), (2.4) 



thus fully including the correlations among the fitted parameters. The uncertainties on 
the nuclides yields were finally obtained with a generalization of the linear propagation 
method described in [|3J[7J|B|. 
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Or conservatively the total ones, if and Oi k were unavailable in a separate form. 
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Figure 1. Comparison between our re- 
gression method (solid line) and the 
one in ref [ H] (dashed line) for the 
2 H(d,p) 3 H reaction. 



Figure 2. As in Fig^ but for the 
2 H(d,n) 3 He process. Different symbols 
are used to denote different data sets. 



3. From the rates to the predictions 

In Fig. El (@J) we plot the net contribution to the right-hand side of the leading and 
main sub-leading reactions to the synthesis and destruction of 7 Li ( 7 Be), obtained for the 
typical value uo^ = 0.023. It is easily seen that the first stage of 7 Li production around 
T ~ 70 keV mainly proceeds through the 4 He(t,7) 7 Li reaction, but 7 Li is soon burned 
via 7 Li(p,a) 4 He reaching a low plateau value. The final 7 Li content is thus mainly given 
by the late 7 Be EC decay: this isotope, synthesized via the key route 4 He( 3 He,7) 7 Be, is 
less easily destroyed through the (experimentally poorly known) channels 7 Be(n,a) 4 He 
and 7 Be(d,p)2 4 He 3 , or the two steps process 7 Be(n,p) 7 Li + 7 Li(p,a) 4 He, now determined 
with good accuracy. The former two reactions are indeed responsible for the bulk of 
the 7 Li error budget, at least if a conservative uncertainty of one order of magnitude is 
assumed. 

We also found that sub-leading processes as 7 Li+p— > 8 Be* — > 7+2 4 He, often neglected 
in BBN studies, for which new measurements exists, are truly marginal, even if their little 
contribution to the final error is similar to the much widely treated 4 He(t,7) 7 Li. The same 
is true e.g. for the other overlooked reactions 6 Li(d,p) 7 Li and 6 Li(d,n) 7 Be. 

We checked that only a handful of reactions dominate the error budget, and, apart for 
useful measurements of poorly known cross sections, a determination of both the magnitude 
and the shape of the 2 H(d,n) 3 He, 2 H(d,p) 3 H, 4 He( 3 He,7) 7 Be reactions at the 1% accuracy 
level over all the interval of interest for the BBN (say, up to ~ 2 MeV) could significantly 



3 The latter has been measured at Louvain la Neuve, see the contribution of C. Angulo to these proceed- 
ings. 
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Leading processes for production and destruction of 7 I_1 



Leading processes for production and destruction of 7 Be 
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Figure 3. The r.h.s. of the Boltzmann 
equations for the reactions relevant for 
the synthesis and destruction of 7 Li. 



Figure 4. As in FigEJ but for 7 Be. Like 
in FigOl the nuclide abundance evolu- 
tion is also shown. 



improve the reliability of the predictions of both 2 H and 7 Li. In these cases, indeed, the 
systematics coming from several experiments dominate the uncertainty, and more reliable 
data are needed, since even very detailed regression methods may fail in these cases. On 
the other hand, the BBN theory would surely benefit of refined studies of the 4 He and 
7 Li observational systematics, as well as from an increase in the statistics of the 2 H/H 
determinations in the high-z damped Ly-a absorbtion systems. 

We also confirm that the 4 He error is dominated by that on the neutron lifetime, while 
the new value of Gn quoted in [E] makes its uncertainty of no relevance for the BBN. 

In summary, we presented some highlights on a new method of data regression and 
a reanalysis and update of the BBN nuclear network. Some differences with the results 
currently quoted in the literature were discussed, but typically a reassuring agreement 
with the usual results has been found, confirming the robustness of the BBN predictions. 

Obviously, the new compilation produced will also turn to be useful to deepen our 
insights on several non standard BBN scenarios. 
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